%作者 : lws
%时间 : 2021/12/19
%%%%good
bs = [0 1 -1];% 分子
as = [1 3 2];% 分母

w = 0:0.01:100;

H = freqs(bs,as,w);

subplot(121);%画图
plot(w,abs(H));
xlabel('\omega'),ylabel('|H(j\omega)|');
title('幅频响应曲线');

subplot(122);%画图
plot(w,angle(H));
xlabel('\omega'),ylabel('\phi(\omega)');
title('相频响应曲线');

%good
t = 0:0.01:10;
[z,p,k] = tf2zp(bs,as);
fvtool(bs,as,'polezero') %画图
legend('零点','极点'),title('极零图')

sys = tf(bs,as);
y = impulse(sys,t);
plot(t,y); %画图
xlabel('t(s)')
title('系统单位冲激响应');
